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{^J , Summary. We study the gravitational clustering of big bang relic neutrinos onto 

. existing cold dark matter and baryonic structures within the flat ylCDM model. 

We then discuss the implications of clustering for scattering-based relic neutrino 
, detection methods, ranging from flux detection via Cavendish-type torsion balances, 

D ■ to target detection using accelerator beams and cosmic rays. 

Q 



1 Introduction 

The standard big bang theory predicts the existence of 10 87 neutrinos per 
flavour in the visible universe. This is an enormous abundance unrivalled by 
1^ \ any other known form of matter, falling second only to the cosmic microwave 

£S| ' background (CMB) photon. Yet, unlike the CMB which boasts its first de- 

\ tection in the 1960s and which has since been observed and its properties 

measured to high accuracy in a series of experiments, the relic neutrino con- 
tinues to be elusive in the laboratory. The chief reason for this is of course the 
feebleness of the weak interaction. The smallness of the neutrino mass also 
makes momentum-transfer-based detection methods highly impractical. At 
present, the only evidence for the relic neutrino comes from inferences from 
other cosmological measurements, such as big bang nucleosynthesis, CMB 
and large scale structure (LSS) data (e.g., [2]). Nevertheless, it is difficult to 
accept that these neutrinos will never be detected in a more direct way. 

In order to design feasible direct, scattering-based detection methods, a 
precise knowledge of the relic neutrino phase space distribution is indispens- 
able. In this connection, it is important note that an oscillation interpretation 
of the atmospheric and solar neutrino data (e.g., [3]) implies that at least two 
of the neutrino mass eigenstates are nonrelativistic today. These neutrinos 
are subject to gravitational clustering on existing cold dark matter (CDM) 
and baryonic structures, possibly causing the local neutrino number density 
to depart from the standard value of n v = hy ~ 56 cm -3 , and the momentum 
distribution to deviate from the relativistic Fermi-Dirac function. 

In this talk, we describe a method that allows us to study the gravitational 
clustering of relic neutrinos onto CDM/baryonic structures. We calculate the 
present day neutrino overdensities in general CDM halos and in the Milky 
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Way, and then discuss their implications for scattering-based relic neutrino 
detection methods — from flux detection via Cavendish-type torsion balances, 
to target detection using accelerator beams and cosmic rays. 



2 Vlasov equation 

The standard procedure for any clustering investigation involving gravity 
only is to solve the Vlasov, or collisionless Boltzmann, equation (e.g., [4, 5]), 

?li= d A +x . d A +i> . d A =Q . (1) 
Dt dr dx dp 

The single-particle phase density fi(x,p, r) is defined so that dNi = fi d 3 x d 3 p 
is the number of i type particles (e.g., CDM, neutrinos) in an infinitesimal 
phase space volume element. The variables x = r/a(t), p = arriiX, and 
dr = dt/a(t) are the comoving distance, its associated conjugate momentum, 
and the conformal time respectively, with a as the scale factor and rrii the 
mass of the ith particle species. All temporal and spatial derivatives are taken 
with respect to comoving coordinates, i.e., '= d/dr, V = d/dx. 2 

In the nonrelativistic, Newtonian limit, equation (1) is equivalent to 

dfi , P dfi dfi 

— + amiVcj) ■ — = 0, (2) 

or arm ox op 

with the Poisson equation 

V 2 <P = ^Ga 2 J2p i (r)6 i (x,r), (3) 



6i(x,r) = P% ( X,T } - 1, pi(x,T) = ^ [ d 3 p fi(x,p,r), 
Pi\ T ) a J 



(4) 



relating the peculiar gravitational potential 4>(x, r) to the density fluctuations 
5i(x,r) with respect to the physical mean pi(r). 

The Vlasov equation expresses conservation of phase space density fi 
along each characteristic {x(t),p(t)} given by 

dx p dp _ , . . 

— = — , -f = -amiV<f>. (5) 
dr ami dr 

The complete set of characteristics coming through every point in phase space 
is thus exactly equivalent to equation (1). It is generally not possible to follow 
the whole set of characteristics, but the evolution of the system can still be 
traced, to some extent, if we follow a sufficiently large but still manageable 
sample selected from the initial phase space distribution. This forms the basis 
of particle-based solution methods. 



2 Unless otherwise indicated, comoving spatial and temporal quantities are used 
throughout the present work. Masses and densities, however, are always physical. 
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3 Solution method and halo density profiles 

A "first principles" approach to neutrino clustering requires the simultane- 
ous solution of the Vlasov equation (1) for both CDM and neutrinos. This 
is usually done by means of multi-component TV-body simulations. In our 
treatment, however, we make two simplifying approximations: 

1. We assume only the CDM component p m contributes to <p m the Pois- 
son equation (3), and p m to be completely specified by halo density 
profiles from high resolution /1CDM simulations [6]. The neutrino com- 
ponent is treated as a small perturbation whose clustering depends on 
the CDM halo profile, but is too small to affect it in return. This as- 
sumption is well justified, since, on cosmological scales, LSS data require 
Pv/Pm = ftv/Om < 0.2 [2]. On cluster/galactic scales, neutrino free- 
streaming ensures that p v j ' p m always remains smaller than Q v jfl m [7]. 

2. Given that assumption f. holds, it follows that not only will the CDM 
halo be gravitationally blind to the neutrinos, the neutrinos themselves 
will also have negligible gravitational interaction with each other. 

These approximations together allow us to track the neutrinos one at a time 
in N independent simulations, instead of following N particles simultaneously 
in one single run. We shall call this "iV-one-body simulation" [1]. 

For the halo density profiles, we use the "universal profile" advocated by 
Navarro, Frenk and White (hereafter, NFW) [8, 9], 

Phaio(r-) = , , u Ps — r -^- (6) 
(r/r s )(l + r/r s ) z 

The parameters r s and p s are determined by the halo's virial mass M v j r and 
a dimensionless concentration parameter c = r v - lr /r s , where r vlr is the virial 
radius, within which lies M v ; r of matter with an average density equal to <5th 
times the mean matter density p m at that rcdshift, i.e., 

4?r. _ o o 4-7T _ o 
M vil = — T HP m a°r* iT = —5 T B.p m fiK h , (7) 

where p m .o is the present day mean matter density. The factor 5th is the 
overdensity predicted by the dissipationless spherical top-hat collapse model, 

18tt 2 + 82y - 39y 2 



y = n m (z) - 1, (8) 



Q m {z) 

with fi m (z) = f2 m ,o/(M m ,o + ^A,oa 3 ) [10]. 

Furthermore, halo concentration correlates with its mass. At z = 0, the 
trend 

/ m . \-°- 13 

was found in [11]. In addition, for a fixed virial mass, the median concen- 
tration parameter exhibits a redshift dependence of c(z) ~ c(z — 0)/(l + z) 
between z = and z = 4. 
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Fig. 1. Relic neutrino number density per flavour, n„ = no, normalised to n v = 
n p ~ 56 cm -3 , for the indicated neutrino and halo virial masses. Results from N- 
one-body simulations are denoted by red (solid) lines. Dotted lines correspond to 
overdensities calculated with the linear approximation. 



4 Clustering in NFW halos 

Using the NFW halo profile (6) as an input, we find solutions to the Vlasov 
equation in the limit p v <C p m - The CDM distribution is modelled as follows: 
We assume a uniform distribution of CDM throughout space, with a spherical 
NFW halo sitting at the origin. For the neutrinos, we take their initial distri- 
bution to be the homogeneous and isotropic Fermi-Dirac distribution with no 
chemical potential. The initial redshift is taken to be z = 3, since, at higher 
redshifts, a sub-eV neutrino has too much thermal velocity to cluster effi- 
ciently. The cosmological parameters used are {Q m , Qa> h} — {0.3, 0.7, 0.7}. 
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Fig. 2. Mass density ratio p v j ' p m normalised to the background mean p„/p m ob- 
tained from iV-one-body simulations for the indicated neutrino and halo masses. 

We solve the Vlasov equation using iV-one-body simulations, as well as a 
semi-analytical linear method. 3 The essential features of the results (Figures 
1 and 2) can be understood in terms of neutrino free-streaming, which causes 
n v /n v to flatten out at small radii, and the mass density ratio p v j ' p m to drop 
substantially below the background mean. Both n v /n v and p v j ' p m approach 
their respective cosmic mean of 1 and p v j 1 p m at large radii. 

Furthermore, we find that the linear method systematically underesti- 
mates the neutrino overdensities over the whole range of halo and neutrino 
masses considered here. Reconciliation with iV-one-body simulations can only 
be achieved if we impose a smoothing scale of > 1 Mpc, or if n v jn v < 3 ~ 4. 



3 The linear approximation [12] consists of replacing df /dp with dfo/dp in (1), 
where /o is the unperturbed Fermi-Dirac function. 
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Fig. 3. Relic neutrino number density per flavour, n v = no, in the Milky Way for 
various neutrino masses. All curves are normalised to n„ = fiy ~ 56 cm~ 3 . The top 
curve in each plot corresponds to the MWnow run, and the bottom to the NFWhalo 
run. The enclosed region represents a possible range of overdensities at z — 0. 

This finding is consistent with the standard lore that perturbative methods 
fail once the perturbations exceed unity and nonlinear effects set in. 

5 Clustering in the Milky Way 

In order to calculate the neutrino overdensity in the Milky Way and, espe- 
cially, their phase space distribution at Earth (r® ~ 8 kpc from the Galactic 
Centre), we need, in principle, to know the complete assembly history of the 
Milky Way. Theory suggests that the galactic bulge and disk grew out of 
an NFW halo via baryonic compression [13, 14]. Our strategy, then, is to 
conduct two series of simulations, one for the present day Milky Way mass 
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Fig. 4. Momentum distribution of relic neutrinos at for various neutrino masses. 
The red (solid) line denotes the MWnow run, while the dashed line represents the 
NFWhalo run. The relativistic Fermi-D irac function is indicated by the dotted line. 
The escape velocity « csc = -\/2|0(r®) is 490 km s" 1 and 450 km s _1 for MWnow 
and NFWhalo respectively, corresponding to "escape momenta" y csc = m u v csc /T^ t o 
of (5.9, 4.4, 3.0, 1.5) and (5.4, 4.1, 2.7, 1.4) for m v = (0.6, 0.45, 0.3, 0.15) eV. 

distribution (MWnow) [15, 16] which we assume to be static, and one for the 
NFW halo (NFWhalo) that would have been there, had baryon compression 
not taken place. The real neutrino overdensity should then lie somewhere be- 
tween these two extremes. Figure 3 shows the possible ranges of overdensities 
at z = 0. 

In all cases, the final momentum distribution at r ffi is almost isotropic, 
with a zero mean radial velocity (v r ), and second velocity moments that 
satisfy approximately the relation 2(vf) = (v^,). Hence, we plot the coarse- 
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grained phase space densities / (Y e , p) only as functions of the absolute ve- 
locity. 

The coarse-grained spectra in Figure 4 show varying degrees of deviation 
from the relativistic Fermi-Dirac function, but share a common feature that 
/ ~ 1/2 up to the momentum state corresponding to the escape velocity from 
the Milky Way at r ffi . This agrees with the requirement that the final coarse- 
grained density must not exceed the maximal initial fine-grained distribution, 
/ < max(/ ) [17, 18, 19, 20, 21]. For neutrinos, max(/ ) = 1/2 at p = 0. 
Thus, our / not only satisfies but completely saturates the bound up to p CS c, 
forming a semi-degenerate state that can only be made denser by filling in 
states abovep csc . 4 



6 Relic neutrino detection 
6.1 Flux detection 

The relic neutrinos' low average momentum (p) = (y) T v $ corresponds to a 
de Broglie wavelength of macroscopic dimension, A" = l/(p) = 0.12 cm/(y). 
Therefore, one may envisage scattering processes in which many target atoms 
act coherently [22, 23] over a macroscopic volume A" 3 , so that the elastic 
scattering rate is proportional to the square of the number of target atoms 
in A -3 . Compared to the case where the neutrinos are elastically scattered 
coherently only on the individual target nuclei, the new rate is enhanced by 
a factor of 

A V A J W cm / \0.1 cm J v ; 

where Na is the Avogadro constant, A is the atomic mass, and pt is the mass 
density of the target material. 5 

Exploiting this effect, a practical detection scheme for the local relic neu- 
trino flux is based on the fact that a test body of density p t at Earth experi- 
ences a neutrino wind force through random scattering events, leading to an 
acceleration given, for Dirac neutrinos, by [22, 23, 27, 28] 

E47T 2 q 
rhyVrei -^N A p t r t o vN 2 m„ v TC i 

flux mom. transfer 

\n v J \ v TC i ) Vg/cnrV \n/{m v v Te \)J 

4 This degeneracy should not be confused with that arising from the Pauli exclu- 
sion principle. 

5 In the case of coherent scattering, it is possible, in principle, to measure also the 
scattering amplitude itself [24, 25, 26], which is linear in Gf- However, a large 
lepton asymmetry is required for a non-negligible effect. 
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with r t < X, <j u n — G 2 F m 2 /ir is the elastic neutrino-nucleon cross section, 
v ici = (\ v ~ v &\) the mean neutrino velocity in the detector's rest frame, 
and u ffi ~ 2.3 x 10 2 km s _1 ~ 7.7 x 10~ 4 c the Earth's velocity through the 
Milky Way. For n v ln v ~ 20, equation (11) gives a t ~ 10~ 26 cm s~ 2 . For 
Majorana neutrinos, a t is further suppressed by a factor of (w TO i/c) 2 — 10 -6 
for an unpolarised target, and v IC \/c ~ 10~ 3 for a polarised one. 

To digest these estimates, we note that the smallest measurable accelera- 
tion at present is > 10~ 13 cm s~ 2 , using conventional Cavendish-type torsion 
balances. Possible improvements with currently available technology to a sen- 
sitivity of > 1CT 23 cm s~ 2 have been proposed [29, 30]. However, this is still 
off the prediction (11) by three orders of magnitude. Therefore, we conclude 
that an observation of this effect will not be possible in the next decade, but 
can still be envisaged in the foreseeable future (thirty to forty years according 
to [32], exploiting advances in nanotechnology) , if our known light neutrinos 
are Dirac particles. Should they turn out, in the meantime, to be Majorana 
particles, flux detection via mechanical forces will be a real challenge. 

Lastly, the background contribution to the acceleration (11) from the solar 
pp neutrinos [flux - 10 11 cm^s" 1 , (E v ) ~ 0.3 MeV (e.g., [31])], < sun ^ 
10~ 27 cm s~ 2 [27], may be rejected by directionality. The background from 
weakly interacting massive particles (WIMPs x, with mass m x ) [27], 

a^ lMP ~ n x v rc \ NaA(j x n 2m x v rc i (12) 

flux mom. transfer 

~ 6 X 10_29 (o.3 GeV/cm 3 ) (wH) ( 100 ) ( lO^cm 2 ) Cm ^ 

should they be the main constituent of galactic dark matter with mass density 
p x = n x m x ~ 0.3 GeV cm~ 3 at r®, can be neglected as soon as the WIMP- 
nucleon cross section <t x n is smaller than ~3x 10~ 45 cm 2 . This should be 
well established by the time relic neutrino direct detection becomes a reality. 



6.2 Target detection 

Detection methods based on the scattering of extremely energetic particles 
(accelerator beams or cosmic rays) off the relic neutrinos as a target take 
advantage of the fact that, for centre-of-mass (cm.) energies, 

V5 = V2m v E hcam * 4.5 (^) (|^) V2 MeV, (13) 

just below the W- and Z-resonances, the weak interaction cross sections grow 
rapidly with the beam energy -Ebcam- 

At accelerators Target detection using accelerator beams does not seem vi- 
able. For a hypothetical beam energy of 10 7 TeV and an accelerator ring of 
ultimate circumference L ~ 4 x 10 4 km around the Earth, the interaction 
rate is roughly one event per year. See [1] for details. 
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With cosmic rays It was pointed out by Weiler [33, 34] (for earlier sugges- 
tions, see [35, 36, 37, 38, 39]) that the resonant annihilation of extremely 
energetic cosmic neutrinos (EECf) — with E > 10 20 eV — with relic anti- 
neutrinos (and vice versa) into Z-bosons appears to be a unique process 
having sensitivity to the relic neutrinos. On resonance, 



the associated cross section is enhanced by several orders of magnitude, 



leading to a "short" mean free path l v = (n v (ff aml )) _1 ~ 1.4 x 10 5 Mpc 
which is only about 48 h times the Hubble distance. Neglecting cosmic evo- 
lution effects, this corresponds to an annihilation probability for EECV from 
cosmological distances on the relic neutrinos of 2 

The signatures of annihilation are (i) absorption dips [33, 34, 40] (see also 
[41, 42, 43]) in the EECV spectrum at the resonant energies, and (ii) emission 
features [44, 45, 46, 47, 48] (Z-bursts) as protons and photons with energies 
spanning a decade or more above the Greisen-Zatsepin-Kuzmin (GZK) cutoff 
at E GZK ~4x 10 19 eV [49, 50]. This is the energy beyond which the CMB 
is absorbing to nucleons due to resonant photopion production. 6 

The possibility to confirm the existence of relic neutrinos within the next 
decade from a measurement of the aforementioned dips in the EECV flux 
was recently investigated in [40] . Presently planned neutrino detectors (Pierre 
Auger Observatory [55], IceCube [56], ANITA [57], EUSO [58], OWL [59], 
and SalSA [60]) operating in the energy regime above 10 21 eV appear to be 
sensitive enough to lead us, within the next decade, into an era of relic neu- 
trino absorption spectroscopy, provided that the EEC^ flux at the resonant 
energies is close to current observational bounds and the neutrino mass is 
> 0.1 eV. In this case, the associated Z-bursts must also be seen as post- 
GZK events at the planned cosmic ray detectors (Auger, EUSO, and OWL). 

What are the implications of relic neutrino clustering for absorption and 
emission spectroscopy? Firstly absorption spectroscopy is predominantly sen- 
sitive to the relic neutrino background at early times, with the depths of the 
absorption dips determined largely by the higher number densities at z >• 1. 
Since neutrinos do not cluster significantly until z < 2, clustering at recent 
times can only show up as secondary dips with such minimal widths in en- 
ergy [61] that they do not seem likely to be resolved by planned observatories. 

On the other hand, emission spectroscopy is directly sensitive to the relic 
neutrino content of the local universe (z < 0.01 ^ tqzk < 50 Mpc). However, 

6 The association of Z-bursts with the mysterious cosmic rays observed above 
Egzk is a controversial possibility [44, 45, 46, 47, 48, 51, 52, 53, 54]. 
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ds ^( s )/ m z - ZkVzGf ^4x 10~ 32 cm 2 



(15) 
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Fig. 5. "Large scale" overdensities (i — v, CDM) in the local universe, with the 
Milky Way at r = 0. The black (solid) line corresponds to the local CDM distribu- 
tion inferred from peculiar velocity measurements [62] (see also [63]) smeared over 
the surface of a sphere with radius r [48] . The dotted line is the neutrino overdensity 
for m v = 0.6 eV, short dash 0.3 eV, long dash 0.15 eV, and dot-dash 0.04 eV. 

since the neutrino density contrasts approximately track those of the under- 
lying CDM above the neutrino free-streaming scale fc^ 1 , it is clear that there 
cannot be a substantial neutrino overdensity over the whole GZK volume 
(~ ^gzk)- Indeed, given the local CDM distribution inferred from peculiar 
velocity measurements (smeared over ~ 5 Mpc), we estimate the correspond- 
ing neutrino overdensity to be < 2 (Figure 5) . Hence the overall emission rate 
cannot be significantly enhanced by gravitational clustering. 

Another possibility is to exploit the fact that there are several galaxy 
clusters (> 10 14 M©) within the GZK zone with significant neutrino clustering. 
One could then search for directional dependences in the emission events 
as a signature of EEC^-relic v annihilation. For example, AGASA has an 
angular resolution of ~ 2° [64] . This is already sufficient to resolve the internal 
structures of, say, the Virgo cluster (distance ~ 15 Mpc, M V ir ~8x 10 14 M Q ) 
which spans some 10° across the sky. From Figure 1, the average neutrino 
overdensity along the line of sight towards and up to Virgo is estimated to 
be ~ 45 and ~ 5 for m v = 0.6 eV and 0.15 eV respectively, given an angular 
resolution of ~ 2°. The corresponding increases in the number of events 
coming from the direction of the Virgo cluster relative to the unclustered 
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case, assuming an isotropic distribution of EECV sources, are given roughly 
by the same numbers, since protons originating from ~ 15 Mpc away arrive 
at Earth approximately unattcnuated. The numbers improve to ~ 55 and 
~ 8 respectively with a finer ~ 1° angular resolution. 

7 Conclusion 

We have conducted a systematic and exhaustive study of the gravitational 
clustering of big bang relic neutrinos onto existing CDM and baryonic struc- 
tures within the flat ylCDM model. Our main computational tools are (i) a 
restricted, iV-one-body method, in which wc neglect the gravitational interac- 
tion between the neutrinos and treat them as test particles moving in an ex- 
ternal potential generated by the CDM/baryonic structures, and (ii) a semi- 
analytical, linear technique, which requires additional assumptions about the 
neutrino phase space distribution. In both cases, the CDM/baryonic gravi- 
tational potentials arc calculated from parametric halo density profiles from 
high resolution TV-body studies and/or from realistic mass distributions re- 
constructed from observational data. 

Using these two techniques, we track the relic neutrinos' accretion onto 
CDM halos ranging from the galaxy to the galaxy cluster variety (M v ; r ~ 
10 12 — ► 10 15 M Q ), and determine the neutrino number densities on scales 
<~ 1 — > 1000 kpc for a range of neutrino masses. We find that the linear 
method systematically underestimates the neutrino overdensities over the 
whole range of halo and neutrino masses considered. Reconciliation with TV- 
one-body simulations can only be achieved if we impose a smoothing scale 
of > 1 Mpc, or if the overdensity is no more than three or four. We there- 
fore conclude that the linear theory does not generally constitute a faithful 
approximation to the Vlasov equation in the study of neutrino clustering on 
galactic and sub-galactic scales (< 50 kpc). However, it may still be useful 
for finding the minimum effects of neutrino clustering in other contexts not 
considered in this work (e.g., the nonlinear matter power spectrum [65]). 

Next we estimate the neutrino phase space distribution in the Milk Way, 
especially in our local neighbourhood at Earth r®, taking also into account 
contributions to the total gravitational potential from the galactic bulge and 
disk. Wc find a maximum overdensity of ~ 20 per neutrino flavour in our 
immediate vicinity, provided that the neutrino mass is at its current up- 
per limit of 0.6 eV. For neutrino masses less than 0.15 eV, the expected 
overdensity from gravitational clustering is less than two. The associated 
coarse-grained momentum spectra show varying degrees of deviation from 
the relativistic Fermi-Dirac function, but share a common feature that they 
are semi-degenerate, with phase space density / ~ 1/2, up to the momentum 
state corresponding to the escape velocity from the Milky Way at r®. This 
means that the neutrino number densities we have calculated here for r ffi 
are already the highest possible, given the neutrino mass, without violating 
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phase space constraints. In order to attain even higher densities, one must 
now appeal to non-standard theories (e.g., [66]). 

In terms of scattering-based detection possibilities, this meager enhance- 
ment in the neutrino number density in the Milky Way from gravitational 
clustering means that relic neutrinos are still far from being detected in fully 
earthbound laboratory experiments. For flux detection methods based on co- 
herent elastic scattering of relic neutrinos off target matter in a terrestrial 
detector, a positive detection could be thirty to forty years away, provided 
that light neutrinos are Dirac particles. For light Majorana neutrinos, another 
~ 10 3 times more sensitivity would be required in the detector for a positive 
signal. Target detection methods using accelerator beams seem equally hope- 
less, unless the accelerator is the size of the Earth and operates at an energy 
of ~ 10 7 TeV. 

Meanwhile, target detection using extremely energetic cosmic neutrinos 
(EECz^, > 10 21 eV) remains the only viable means to confirm the existence of 
big bang relic neutrinos within the next decade or so. Resonant annihilation 
of EECf on relic neutrinos can be revealed as absorption dips in the EECf 
flux (e.g., [40]), or as emission features in the Z-decay products. However, 
since absorption spectroscopy is largely insensitive to late time (z < 2) relic 
neutrino clustering, our findings here have little impact on the conclusions 
of [40]. On the other hand, emission spectroscopy is sensitive to the relic 
neutrino content of the local GZK zone, V G zk ~ 50 3 Mpc 3 . While we find 
no significant large scale clustering within Vqzk and therefore no significant 
enhancement in the overall emission rates, it is still conceivable to exploit 
the considerable neutrino overdensities in nearby galaxy clusters, and search 
for directional dependences in the post-GZK emission events. For the Virgo 
cluster, for example, we estimate the event rate from the central 1° region to 
be ~ 55 and <~ 8 times the unclustered rate for neutrino mass m v = 0.6 eV 
and 0.15 eV respectively, assuming an isotropic distribution of EECV sources. 
Planned observatories such as the Pierre Auger Observatory [55] , EUSO [58] 
and OWL [59] will have sufficient angular resolution to, in principle, see 
this enhancement. However, considering the rapidly improving constraints 
on both the EECis flux and neutrino masses, it remains to be seen if the 
enhancement can indeed be observed with enough statistical significance [67] . 
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